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We study the effect of Majorana transition magnetic moments on the fiavor evolution of neutrinos 
and antineutrinos inside the core of Type-II supernova explosions. We find non-trivial collective 
oscillation effects relating neutrinos and antineutrinos of different flavors, even if one restricts the 
discussion to Majorana transition electromagnetic moment values that are not much larger than 
those expected from standard model interactions and nonzero neutrino Majorana masses. This 
appears to be, to the best of our knowledge, the only potentially observable phenomenon sensitive 
to such small values of Majorana transition magnetic moments. We briefly comment on the effect of 
Dirac transition magnetic moments and on the consequences of our results for future observations 
of the flux of neutrinos of different flavors from a nearby supernova explosion. 

PACS numbers: 14.60.Pq, 97.60.Bw 

I. INTRODUCTION 

The understanding of neutrino flavor oscillations in the presence of a dense neutrino background has improved 
dramatically over the last several years The collective effects mediated by neutrino self-interactions lead to 

new oscillation phenomena, some of which may be imprinted, in an observable way, in the flavor and energy spectra 
of neutrinos from Type-II supernova explosions. 

These novel effects are very sensitive to some of the currently unknown neutrino masses and mixing parameters, 
in particular the neutrino mass hierarchy. For example, in the case of the so-called inverted neutrino mass hierarchy 
[l3l4T5| . it has been established that there is a "flavor spectral split" in the neutrino flux [l6l - [23j . while none is 
observed in the case of the normal neutrino mass hierarchy. Using initial z/g , Pe 7 , i^x neutrino fluxes as a function 
of the neutrino energy from the steep power law model computed in [23 | (with p = 10 and q = 3.5) and depicted in 
Fig. mjleft), and assuming the matter potentials are as depicted in Fig fright), we reproduce these results under a 
two-flavor, single-angle approximation in Fig. [5] See j25l - l32j and [l^, [33l - l35j for three-flavor and multi-angle analyses 
of the problem. Our results, which will be discussed in more detail in Section [TO are in good agreement with previous 
calculations that make use of the same initial fluxes and matter profile [ijj l23l- While the split in the neutrino 
spectrum in the case of an inverted hierarchy is present only for a non-vanishing mixing angle 9, its location and 
shape are not very sensitive to the actual value of 0. This is a testament to the fact that due to the high degree of 
non-linearity there are "switch-on" effects for certain parameters in the Hamiltonian. 

The recent discovery that the neutrino mixing angle ^13 [s^] is not zero j38l - l42| implies that something similar to 
the phenomenon depicted in Fig. [5] is expected from supernova neutrinos. It should be noted, however, that while 
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FIG. 1: (Left) Initial Ve,i>s,i'x,Vx neutrino fluxes as a function of the neutrino energy, and (right) the strength of the matter 
potential and the effective neutrino self-interaction as a function of the distance from the center of the supernova explosion. 
The matter and self-interaction potentials are deflned as v^GfJ^b and \f2GFrLv, respectively. See Section [II] for details. 
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FIG. 2: Initial (dashed, r — 50 km) and final (solid, r — 200 km) i'e,i'e,Vx,i'x neutrino fluxes as a function of the neutrino 
energy for the normal (left) and inverted mass hierarchy (right), for | Am^| = 2.5 x 10"'^ eV^ and sin^ d = 10~^. See [s^l for an 
animation of the radial evolution of the fluxes. 



the energy spectra of neutrinos are very different for the two hierarchies, it is not yet entirely clear whether a robust 
neutrino mass hierarchy measurement with supernova neutrinos is practical considering the uncertainties in the initial 
fluxes [43| . Here we will not concern ourselves with this very important practical issue. 

In this paper we study the effect of nonzero neutrino magnetic moments in the presence of a strong magnetic field 
on collective neutrino oscillations. The effect of neutrino magnetic moments in neutrino oscillations in variable matter 
density is an interesting phenomenon that was first identified as a potential solution to the solar neutrino puzzle 
several years ago [13,143], but its consequences for collective effects in supernova neutrino oscillations are, to the best 
of our knowledge, yet to be discussed in the literature. 

We assume that neutrinos are Majorana fermions. When the standard model is augmented to include nonzero 
Majorana neutrino masses and lepton mixing, electroweak interactions lead to very small but nonzero transition 
magnetic moments. These results are summarized in Section Hill along with our assumptions regarding the magnetic 
field inside the supernova explosion. We limit our discussion to Majorana neutrinos for several reasons. Some 
are related to technical issues concerning numerical calculations; Majorana neutrinos only have nonzero transitional 
magnetic moments and hence there are fewer free parameters, and only the transverse component of the magnetic field 
affects their fiavor evolution. More important, in the case of Majorana neutrinos, the magnetic moment interactions 
with the magnetic field mediate "ncutrino-antineutrino oscillations" and one naively anticipates the possibility of 
spectral splits between neutrinos and antineutrinos of different flavors. Our results, presented in Section ITVl confirm 
this suspicion and reveal that even very small values the magnetic moment, perhaps those close to the "standard 
model" level, may play a significant role in the evolution of neutrino fluxes due to the high degree of non-linearity in 
the problem. 

In our attempt to draw attention to the importance of the transition magnetic moment (and in order to simplify the 
discussion), throughout, we are limiting our discussion to the so-called single angle approximation, and pretend there 
are only two neutrino/ antincutrino flavors: j/g and Vx, the latter a linear combination of and Vr- We comment on the 
inclusion of multi-angle and three- fiavor effects in the conclusions. Section |Vl There we also address other simplifying 
assumptions, and comment on our expectations regarding the impact of magnetic moments if the neutrinos are Dirac 
fermions. 



II. COLLECTIVE NEUTRINO OSCILLATIONS 



We discuss neutrino oscillations inside the supernova using the density matrix formulation. The evolution of the 
neutrino states is given by 



«P = [P,H]. 



(1) 
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We take p to be a 4 x 4 matrix^ which includes both "neutrino" (left-handed) and "antineutrino" (right-handed) 
states, and H is the Hamiltonian, which governs the time evolution of the system. This will prove useful when we 
consider the effect of the transition magnetic moment. 

We define the density matrix p in the flavor basis such that its initial state is diagonal: 



f Pee Pex Pee 

Pxe Pxx Pxe 

Pee Pex Pee 

\Pxe Pxx Pxe 



(2) 



Pij 



are defined as follows: 



(27r)^<S(^)(p-p')p.,(p) = < 



'(a](p)ai(p')) ifi,j = e,x, 

(6](p)6i(p')) ifi,j = e,x, 

(a](p)5i(p')) li i ~ e,x and j = e, x, 

_ (6](p)ai(p')) if i = e,a; and j = e,x. 



(3) 



We use the notation of [i^ and denote the creation and annihilation operators for neutrinos (anti-neutrinos) by a^{b^) 
and a{b) respectively.^ We also define p'^ to be equal to p with a replaced by b and vice-versa, 



P = 



Pee Pex Pee Pex 

Pxe Pxx Pxe Pxx 

Pee Pex Pee Pex 

\Pxe Pxx Pxe Pxx/ 



The density matrix can be obtained from p by replacing the blocks 11 22 and 12 o 21. 

The Hamiltonian H will be divided into three components: H = Hyac + f^mat + Hseif- Hseif^ 
neutrino-neutrino self-interactions, is 



(4) 



which contains the 



H. 



self 



V2G 



dE GHp{E) p{EY*)G + ictTr {{p{E) p{Ey*)G) , 



(5) 



where n,^ is the total effective number of neutrinos (number of neutrinos plus the number of antincutrinos) , G^? is the 
Fermi constant, and G is a matrix of dimensionless coupling constants. 



G = 



/I 
10 
0-10 

Vo -I, 



(6) 



A derivation of Eq. [5l along with some of the notation, is discussed in Appendix [X] Note that we have a complex- 
conjugation in the second component of Eq. ([5]),'^ which is different from the convention of p defined in (46j . 
The interactions of neutrinos with the rest of the environment are contained in 



/^/2GF(r^e 



V 



1 I 





-V2Gf- 








V2GF(ne 




2 / 



(7) 



\/2Gf^I 



This term is the same as the usual MSW term, written in the form of a 4 x 4 matrix. ne,„ are, respectively, the 
electron and neutron number densities. The neutral current interactions can in principle play a role in the presence 



^ In the absence of effects that mix neutrinos and antincutrinos, it is customary (see, for example, [Zsll ). and more convenient, to describe 
two-flavor oscillations, including antincutrinos, with a 2 X 2 density matrix. 

More precisely, we also extend the notation in [4(| to include the "neutrino-antineutrino" elements (block-off-diagonal components) of 
the density matrix. 

^ Using the fact that p is Hermitian, it is easy to see that p*g = pee, pja, = Pa;e, etc. 
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of Icpton number violation as there is the possibihty of neutrino-antincutrino oscillations. Henceforth, however, we 
neglect them as they do not qualitatively impact our results and would only render the discussion less clear. 

Finally, the relation between the initial states and the vacuum Hamiltonian eigenstates is contained in the "vacuum" 
part of the Hamiltonian, given by (in the 4x4 form used here), in the absence of magnetic fields. 



/-w cos 26* a; sin 26* 
a;sin26' a;cos26' 

cos 26* w sin 26* 

V cjsin26l wcos26'y 



(8) 



where w = is the vacuum oscillation frequency, and is the mixing angle that relates the weak and flavor 

bases. For all practical purposes, 9 is to be interpreted as the lepton mixing angle 6*13 and Am^ as Amfg ^ Amjg 
371 ]. Throughout, unless otherwise noted, we fix 6* = 10^'^ radians, in spite of the recent evidence for a much larger 
^13. This is done in order to focus on the switch-on effect due to collective oscillations and in order to suppress any 
additional structure due to a "large" mixing angle. Figure [5] was generated by numerically solving Eq. ([1]) for the 
4x4 density matrix as a function of time using the initial conditions and matter profiles depicted in Figure [T] 

In the presence of a nonzero Majorana transition magnetic moment and a background magnetic field, the "vacuum" 
Hamiltonian is modified (see, e.g., [47j): 



'UJ cos 29 UJ sin 29 fiB- 



T 

_ I u]sin29 a;cos26' —fiBr . ..^ 
-fiBx cos 26* w sin 261 I ' 
IiBt UJ sin 29 uj cos 29/ 

where Bt is the component of the magnetic field transverse to the neutrino momentum, and /i is the magnitude of 
the Majorana neutrino magnetic moment. The Majorana transition magnetic moment matrix is antisymmetric and, 
in the case of two neutrino flavors, it is parameterized by only one number. 



III. NEUTRINO MAGNETIC MOMENTS AND MAGNETIC FIELDS INSIDE SUPERNOVAE 

When neutrinos are massive, weak interactions lead, at one-loop, to nonzero neutrino magnetic moments. Assuming 
only standard model interactions, if there were only one massive Dirac neutrino, its magnetic moment would be (for 
a recent overview, see [47| ) 

ieGprriy ( ( mj 



87r2^/2 ' ^ V^^^. 
= 3.2 X 10-1^ (^) MB, (10) 

where mg is a charged lepton mass, Mw is the ly-boson mass, and iib = 5.788 x 10~^ eV/gauss is the Bohr magneton. 
For a neutrino mass of order a tenth of an eV, fin ^ 3 x lO^^^/xs. For more than one neutrino species, the transition 
magnetic moments are GIM suppressed with respect to the diagonal ones. Taking into account that all mixing angles 
in the lepton mixing matrix are large, the GIM suppression factor is of order ni^/M^ '--^ 5 x 10""*. 

For Majorana neutrinos only transition magnetic moments are allowed: the magnetic moment matrix is totally 
antisymmetric, fifj ~ ^l^ji ~ 1,2,3. The magnitudes of the Majorana neutrino transition electromagnetic 

moments depend on the neutrino Majorana phases (or, if CT'-invariancc is conserved, on the relative CT^-paritics of the 
different neutrino mass eigenstates). For generic Majorana phases one can estimate that standard model interactions 
lead to, assuming all neutrino masses have the same order of magnitude = 0.1 eV, /dfj = ^l^ji ^ lO^^/^z? ^ 
10~2*/iB, where fio is given by Eq. (fTO|) . For more details see, for example, [l^. In theories beyond the standard 
model, much larger values for the Majorana neutrino transition magnetic moments are, of course, allowed. Our current 
understanding of neutrinos interactions constrains /i < lO^^^^s [37[, some nine orders of magnitude larger than fijj 
above. 

The Hamiltonian Eq. ^ depends on the combination ^Bt and hence we need to discuss the expected magnetic 
fields inside the supernova explosion (we drop the transverse subscript T on B hereafter) . We make use of a toy model 
for the magnetic field and assume it falls off like l/r^, where r is the distance from the center of the supernova, and 
assume a magnitude of lO"'^^ gauss at 50 km, which is consistent with estimates in the literature (see, for example, 
[isl p and might be considered, by some authors, to be conservative (see, for example, [ioj): 

„/ N IT /50 km\ ^ 

B{r) = 10^2 gauss. (11) 
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It should be noted that the magnetic field in the supernova is probably toroidal in nature. Depending on the direction 
of neutrino propagation, it may decrease more gradually than assumed here, or even increase in magnitude at small 
radii [50l |. 

One needs also address the orientation of the magnetic field with respect to the supernova-earth direction, since it 
is not just the magnitude of the magnetic field but also its direction that affects the evolution of the neutrino fiuxes. 
The magnetic field has no effect if it is parallel to the direction of neutrino propagation. Furthermore, we do not 
consider twisting magnetic fields. A dipolar magnetic field is sufficient to illustrate the importance of the transition 
magnetic moment in collective oscillations. Finally, it is also overly simplistic to consider a magnetic field which is 
fairly smooth without any turbulence, but the addition of turbulence can only make things more interesting. 

It will prove useful to define 

if^DB)sMir)^^,nBo[^-^j ^ 1.9 x IQ-^" (^— — j (^^j , (12) 

for Bq = 10^^ gauss and = 3.2 x 10^^'^ piB- We are interested in positions 200 km > r > 50 km. Our results 
will be presented in units of (^jLdB)sm, including the dependency. In these units, standard-model-like Majorana 
neutrino magnetic moments, assuming our estimates for the magnetic field are appropriate, correspond to fJ.B{r) ~ 
10~^(/i£)i?)sM- Magnetic field orientation issues may mean smaller effective ^B values, while magnetic fields larger 
than the ones considered here mean larger fiB. If there is new physics beyond the standard model, much larger values 
— many orders of magnitude larger! — are allowed. 



IV. RESULTS 



When the neutrino magnetic moment effects are negligible, collective neutrino oscillations, for the values of the 
parameters discussed in Sec. |lll are depicted in Fig. These results are well understood in the literature and can be 
qualitatively described, for the inverted mass hierarchy, as a "spectral swap" between i/g and for energies between 
(roughly) 5 MeV and 18 MeV. 

Technically, the addition of transition magnetic moment significantly increases the complexity of the differential 
equations, to the extent that it is very challenging to simply "code" them without the possibility of typographical 
errors. To circumvent this, we use SAGE [5l| to automatically generate the C-code for the input functions. We per- 
formed multiple checks in order to verify that our results are, numerically, reliable, including exploring circumstances 
when either fiB or 9 = vanish. The former check has already been discussed. Eq. (jH]) reveals that a similar behavior 
is to be expected in the case fiB ^0,9 = 0, except that, this time, the "swapping states" are z/g and (and their 
antiparticles).* Our results for 9 = and ^B{r) = 1Q~^{^dB)sm are depicted in Figure |3] and clearly show that 
the expectations are confirmed, this time for a normal neutrino mass hierarchy. Like in the standard case {jiB = 0, 
9 0), one can verify the conservation of 'lepton number' in the sense that Ai/g is equal to Ai?e. Indeed, this is one 
of the "sanity checks" we performed in order to establish the reliability of the numerical results presented here. 



* In the upcoming figures, this is not immediately obvious due to the fact that and Ux have identical initial spectra. 
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FIG. 4: Initial (dashed, r — 50 km) and final (solid, r — 200 km) i'e,Vs,i'x,i'x neutrino fluxes as a function of the neutrino 
energy, including the effect of transition magnetic moment for neutrinos (left) and antineutrinos (right), for a normal mass 
hierarchy and i-iB{r) = 10^^ {i-IdB)sm (see Sec lIII|) . See [s^ for an animation of the radial evolution of the fluxes. 
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FIG. 3: Initial (dashed, r — 50 km) and final (solid, r = 200 km) VejVsjVxjVx neutrino fluxes as a function of the neutrino 
energy for 9 — and fiB{r) — 10^^(/^d-B)sm (see Sec lIII[) . for a normal (left) and inverted (right) neutrino mass hierarchy. 
See [391 for an animation of the radial evolution of the fluxes. 



We draw attention to the fact that ^B{r) = 10 ^{fJ.DB)sM is larger (by a couple of orders of magnitude) than 
the "baseline" standard model expectation for massive Majorana neutrinos discussed in Sec. IIIII Numerically, we find 
that /iB(r) = IQ-'^i^ioB) SM does not lead, given the initial fluxes considered here, to a significant effect. We return 
to this issue in the next section. 

Naively, one could anticipate that when both 9 and fiB arc nonzero, the combined collective oscillation effects 
would "add up," and the final fiux could be obtained by starting with the standard case {^B = 0) and applying the 
logic mentioned earlier in order to estimate the effect of the nonzero ^B. Our results for sin^ 6* = 10""*, nB{r) = 
10~^(/iDi?)sM are depicted, for a normal and inverted neutrino mass hierarchy, in Figs.|4l[5l The figures reveal more 
structure in the final neutrino fluxes than naively anticipated. There are effects for both the normal and inverted mass 
hierarchies, and all flavor neutrinos and antineutrinos are affected. As we mentioned earlier, in the absence of the 
transition magnetic moment, the lepton number is conserved at each and every energy. Once the transition magnetic 
moment is "turned on," lepton number is violated at individual energies. Total lepton number (integrated over all 
neutrino energies), however, is conserved. 

Different initial fluxes lead to quantitatively different results. We consider initial fluxes, depicted in Figurejpt^lcft) 
(see [l^) where the ratios L^^ : L^^ and Lp^ : L^^ are signiflcantly higher than the ones in computed in [SJ and 
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FIG. 6: Initial (left, r = 50 km) and final (right, r = 200 km) i'e,i'e,Vx, Vx neutrino fluxes as a function of the neutrino energy, 
assuming different initial neutrino spectra from the previous figures: higher L^^ : L^^ ratio and Lp^ : L^^ ratio, for — (no 
magnetic moment effect) and an inverted neutrino mass hierarchy. For these initial fluxes, collective oscillations are virtually 
absent for the normal neutrino mass hierarchy. See [s^ for an animation of the radial evolution of the fluxes. 



repeat the computations performed above. These results, depicted in Figure IHJright)^ and Figure[71 are quantitatively 
different from the ones in Figs. [U UJ and [5] Here, we see that the nonzero magnetic moment effects are much more 
pronounced — there are multiple splits — in the case of a normal mass hierarchy, while in the case of an inverted one 
the effect of the nonzero magnetic moment is virtually absent. 

The comparison of Figs. E^right) , [7]with Figs. [21 IH and [5] illustrates the challenge of disentangling "particle physics" 
effects from "astrophysics effects." For the two sets of figures the neutrino parameters arc, pairwise, the same, while 
the initial neutrino fluxes arc different. The structure in the final fluxes is clearly different. Regardless of the initial 
neutrino spectra, however, our results allow one to conclude that modest values for the neutrino magnetic moment, in 
the region of the supernova explosion where collective effects are dominant (r < 200 km) , can significantly modify the 
energy spectrum of neutrinos of different flavors, and that, once one combines the "standard" ^13 effect with a nonzero 
magnetic moment, spectral splits — sometimes multiple splits — are expected for both neutrino mass hierarchies. 

In order to relate theoretical fluxes with potential obscrvables at terrestrial experiments, it is necessary to further 
"evolve" the neutrinos from r = 200 km to r — >■ c». Unfortunately, the standard matter effects, which dominate outside 
the "collective oscillations" region, cannot be treated in a satisfactory manner in the two-flavor approximation, as the 
multiple level crossings play an important role at larger radii. This by itself is a very interesting topic [Tol. [52j. and 



^ For this flux and a normal hierarchy, spectra are virtually unaffected by collective oscillations. 
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has to be combined with the study of collective oscillations in order to draw any conclusions regarding the observable 
signatures due to the transition magnetic moment. The three flavor analysis of this problem, which is beyond the 
intentions of this manuscript, would take us a step closer towards possible observable signatures in the final flux. It 
contains additional intricacies, including potentially three distinct transition magnetic moments. A detailed analysis 
of the three flavor effect in collective oscillation, and its interplay with the standard MSW effect in supernovae at 
larger radii, is work in progress which we hope to report in the near future. 



V. CONCLUSION 



We have studied the effect of nonzero Majorana transition magnetic moments on neutrino collective oscillations in 
the interior of supernova explosions. Our main results are illustrated in Figures [3l |4l and [5l and in [36| . We find 
that, even for modest values of the neutrino magnetic moment and the magnetic field inside the supernova, nontrivial 
effects are expected. These are similar in character and magnitude to effects of nonzero leptonic mixing angles, which 
have been thoroughly studied in the literature. The interplay between "magnetic moment effects" and "flavor mixing 
effects" is also interesting and may lead to the presence of nontrivial structure in the neutrino flavor spectra for both 
the normal and the inverted neutrino mass hierarchies. 

We find that nontrivial magnetic moment effects occur for fi values that arc only a couple of orders of magnitude 
larger than the standard model expectations, assuming no new physics beyond light neutrino Majorana masses. It 
is, arguably, virtually impossible to imagine another potential observable — astrophysical or terrestrial — that is, in 
practice, sensitive to such tiny neutrino magnetic moments. Our result is dependent on the very simple, albeit mostly 
conservative, assumptions we made regarding the magnitude, form, and radial dependency of the magnetic field in 
the supernova core. The fact that one is sensitive to unpreccdentcdly small values of the neutrino magnetic moment 
is, however, robust under changes of the hypothesis regarding magnetic fields. 

We have restricted our discussion to the propagation of neutrinos inside the region where collective oscillations 
dominate (r < 200 km). Supernova neutrino spectra, when observed in terrestrial experiments, will be further 
"processed" by neutrino propagation further downstream within the supernova, including regions where the "standard" 
matter effects are strong, and potential propagation through the Earth. Throughout, if the neutrino magnetic moments 
are much larger than standard model expectations (and much larger the the values considered here) other magnetic 
moment effects are expected outside the region under consideration here. Disentangling all of these effects from the 
magnetic-moment induced collective effects discussed here is a formidable task we have not addressed here at all. We 
do, however, point out that the magnetic moment effects considered here may partially mimic effects characteristic of 
the wrong neutrino mass hierarchy, and that the entire picture is also sensitive to the initial neutrino flavor spectra 
inside the supernova core. 

Disentangling the physics contained in measured supernova neutrino spectra (which will happen in the future, 
hopefully at multiple detectors) is a well known challenge. We point out that neutrino electromagnetic properties are, 
potentially, an important component of this "entanglement." The physics opportunities, on the other hand, are very 
exciting. For example, if astrophysical uncertainties are under control, and if the neutrino mass hierarchy is known, it 
is possible that, by observing supernova neutrinos, we may be able to determine that the neutrino magnetic moments 
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are nonzero. Furthermore, if Majorana magnetic moments are qualitatively different from Dirac ones (c/. discussion 
in the last paragraph), it is possible we may be able to indirectly infer that neutrinos are Majorana fermions just by 
observing effects of neutrino magnetic moments from supernovae. Such an inference would be completely independent 
from other observables sensitive to the nature of the neutrino. 

Our main goal was to point out that even very small neutrino Majorana transition magnetic moments can sig- 
nificantly impact the flavor structure of supernova neutrinos. Along the way, we have made several simplifying 
assumptions that need to be revisited if one is to obtain more robust quantitative results. In more detail, we assume 
the existence of only two neutrino flavors and have restricted our analysis to the single angle approximation. There 
may be other interesting features of collective oscillations we have missed due to these approximations. In the ab- 
sence of neutrino magnetic moments, for example, the inclusion of three-flavor effects (see [25l - [3^ ) and multi-angle 
analyses (see [231 [33l - l35| ) have revealed some interesting effects not properly captured by the two- flavor, single-angle 
approximation. 

Throughout, we have restricted our analysis to Majorana neutrinos. If the neutrinos are Dirac fermions the problem 
is both technically more involved and qualitatively different. If the neutrinos are Dirac fermions, the magnetic moment 
interacting with the local magnetic field leads to "active-sterile" oscillations (where the new right-handed neutrino 
degrees of freedom act, for all practical purposes, as sterile neutrinos) and the "two-flavor" approximation is, in reality, 
a "four-flavor" problem. Furthermore, Dirac neutrinos can also carry a diagonal magnetic moment and the number 
of free parameters is, in the two- flavor-case, three (as opposed to one in the Majorana case). Qualitatively, using 
previous results as guidance, we don't expect any sizable effects. The reason is simple. Collective oscillations obey 
certain conservation laws and the fact that the initial flux of sterile neutrinos is zero implies that no "swaps" can take 
place — there are no sterile neutrinos around to trade places with the active ones! Preliminary results indicate that 
this naive picture is indeed realized. 



Appendix A: The Self-interaction Hamiltonian 



At the quantum field level, the neutrino-neutrino interaction Hamiltonian density is 



where /, /' = e,x and 



Hseif = V2GFY,{¥r^^){^^'l>.^^'). (Al) 
/,/' 



V'^(^) = / (a^«(p)e-^^^ + 6^t„(p)e^P-) . (A2) 



The creation and annihilation operators were defined in Eq. Q. The related field ip-^{x) can be defined as Cip^{x)C. 
While we are dealing with Majorana particles, we can still safely consider fields made up of different creation and 
annihilation operators corresponding, respectively, to left-hand and right-handed particles. 

For sake of presentation, we derive a particular element of Hgeif,^ namely {Hseif)ex which, (using d?q = 
d^q/(27r)'^2i?g) can be written as. 



= -V2Gf y dVrfVa;(i/5(p')|ViW'1i'e(p)) 

X (q|^^7^^^|q). (A3) 
Here, the one-particle neutrino states of well-defined momentum are defined as 

|;./(p)) = y2^a^t|o), (A4) 
for / = e, X. These eigenstates are normalised according to 

(p7'Ip, /) = 2E^{2i,f5^^\v - v')5^f'. (A5) 



All other ones can be derived in exactly the same way. 



10 



The same expressions with a b apply for antineutrino states. |q) are a complete set of states (including neutrinos 
and antineutrinos, and all flavors). In Eq. (jA3p we have used the Fierz identity in the second line, so that we can see 
that for the second matrix element the only creation and annihilation operators that survive arc beol and alb^ with 
a relative minus sign between them. Using the elements of the density matrix, 

(q|Vi'7''V'''| q) a pes{i'e\JLa%\iys) + Pxe{'^x\J^^beal\iys), (A6) 

A relative negative sign between the two terms arises from Jj'g = uj^^v ~ ^Jgxi these arc CV conjugates of 
each other. Finally, we use the standard technique of summing over the spins by taking the trace which leads to a 
dependency of (Hseif) on the angle between the momenta of the neutrinos, 

{H,eif)ex = V2Gf j d^q{p,,-p,^)(^^^y (A7) 

= V2Gf J Sq,{p^s- Pex){l-cos{9pq)) . (A8) 

The potential derived above can be averaged over the angle between relative velocities Opq to obtain an average 
potential per particle in the thermal bath. The potential experienced by a particle traveling through the thermal bath 
is proportional to the number of particles in the bath. The above equation can be written in the form presented as 
Eq. [5] with the help of p^: pxe = pTx- 

The matrix elements of (Hseif)ex can be obtained from the expression for {Hseif)ex in [46j as follows. Every 
time we replace a neutrino by its antiparticle, the corresponding current (see Eq. (jA3[) ) acquires a relative negative 
sign. This means that, in order to derive the full self- interaction Hamiltonian, the number of matrix elements that 
need to be calculated is very small. This sign flip is also the reason why the coupling constant matrix has the form 
G ~ diag(l, 1,-1,-1), in Eq. ([5]). This form of the coupling constant matrix ensures that there is a negative sign 
whenever any two flavors in the equations of motion arc replaced by their charge conjugate. 
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